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been  constructed  and  the  improvement  is  documented  by  Measurements  before  and 
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PREFACE 


This  effort  was  conducted  by  Syracuse  University  under  the  sponsorship 
of  the  Rome  Air  Development  Center  Post-Doctoral  Program  for  Rome  Air 
Development  Center.  Mr.  Charles  Pankiewicz,  RADC/RBTV,  was  the  task  project 
engineer  and  provided  overall  technical  direction  and  guidance. 

The  RADC  Post-Doctoral  Program  is  a cooperative  venture  between  RADC 
and  some  sixty-five  universities  eligible  to  participate  in  the  program. 
Syracuse  University  (Department  of  Electrical  Engineering) , Purdue  Univer- 
sity (School  of  Electrical  Engineering),  Georgia  Institute  of  Technology 
(School  of  Electrical  Engineering),  and  State  University  of  New  York  at 
Buffalo  (Department  of  Electrical  Engineering)  act  as  prime  contractor 
schools  with  other  schools  participating  via  sub-contracts  with  prime 
schools.  The  U.S.  Air  Force  Academy  (Department  of  Electrical  Engineering), 
Air  Force  Institute  of  Technology  (Department  of  Electrical  Engineering), 
and  the  Naval  Post  Graduate  School  (Department  of  Electrical  Engineering) 
also  participate  in  the  program. 

The  Post-Doctoral  Program  provides  an  opportunity  for  faculty  at 
participating  universities  to  spend  up  to  one  year  full  time  on  explora- 
tory development  and  problem-solving  efforts  with  the  post -doctoral s 
splitting  their  time  between  the  customer  location  and  their  educational 
institutions.  The  program  is  totally  customer-funded  with  current  pro- 
jects being  undertaken  for  Rome  Air  Development  Center  (RADC),  Space  and 
Missile  Systems  Organization  (SAMSO),  Aeronautical  System  Division  (ASD) , 
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Electronics  Systems  Division  (ESD),  Air  Force  Avionics  Laboratory  (AFAL) , 
Foreign  Technology  Division  (FTD) , Air  Force  Weapons  Laboratory  (AFWL) , 
Armament  Development  and  Test  Center  (ADTC),  Air  Force  Communications 
Service  (AFCS),  Aerospace  Defense  Command  (ADC),  HQ  USAF,  Defense  Com- 
munications Agency  (DCA) , Navy,  Army,  Aerospace  Medical  Division  (AMD), 
and  Federal  Aviation  Administration  (FAA) . 

Further  information  about  the  RADC  Post-Doctoral  Program  can  be 
obtained  from  Mr.  Jacob  Scherer,  RADC/RBC,  Griffiss  AFB,  NY,  13441, 
telephone  Autovon  587-2543,  Commercial  (315)  330-2543. 
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1.  INTRODUCTION 


1.1  Background 

The  Newport  Antenna  Test  Range  of  Rome  Air  Development  Center  is  con- 
cerned with  the  testing  and  evaluation  of  ground-based  and  aircraft-mounted 
antennas. 

At  present,  its  principal  test  sites  are  located  on  Tanner  Hill  and 
Irish  Hill  as  indicated  in  Figure  1. 

A number  of  studies  have  been  conducted  in  the  past  with  the  objec- 
tives of  evaluating  the  performance  of  the  Test  Range  and  recommending  ways 
of  enhancing  its  facilities.  Titles  of  reports  giving  some  of  the  results 
of  these  prior  studies  are  listed  in  Appendix  I. 


1.2  Objectives 

The  major  objectives  of  the  present  study  were  two-fold: 

1. )  Determine  the  location  of  major  reflection  zones  for  the 

existing  short  and  long  test  ranges  and  the  effect  of  reflections  on 

antenna  test  measurements  over  the  frequency  range  from  200  MHz  to  2.0  GHz. 
Recommend  techniques  that  can  reduce  the  reflections  and  analyze  the  ex- 
pected improvement  in  measurement  accuracy  obtained  by  each. 

2. )  Identify  and  analyze  the  site  and  instrumentation  limi- 
tations to  be  anticipated  in  extending  the  Newport  Test  Range  antenna  measure- 
ment capability  from  18  to  AO  GHz. 

The  report  includes  some  well-known  results  which  are  found  in 
standard  references,  but  these  have  been  included  here  for  completeness. 


i 


i 


! 


f 


Figure  1.  Topographic  View  of  the  RADC  Antenna  Test  Site 


The  experimental  results  and  data  Included  in  this  report  were  ob- 
tained by  personnel  of  Rome  Research  Corporation  (RRC)  under  the  supervision 
of  Mr.  Eugene  Nordell.  This  assistance  is  gratefully  acknowledged. 


REFLECTIONS 


2.1  Statement  of  the  Problem 

A specific  antenna  pattern  is  usually  obtained  by  employing  the  antenna 
as  a receiving  antenna.  As  a consequence  of  the  Reciprocity  Theorem1,  the 
transmitting  and  receiving  patterns  of  an  antenna  are  the  same  in  a linear, 
passive,  isotropic  medium,  if  the  illuminating  field  is  a uniform  plane  wave. 
Major  sources  of  disturbances  in  an  otherwise  uniform  field  frequently  are 
reflections. 

Any  ponderable  body  such  as  the  terrain,  man-made  structures,  vegeta- 
tion, etc.  may  give  rise  to  reflections  which  disturb  the  field  in  the  neigh- 
borhood of  the  test  antenna. 


2.2  Factors  Affecting  Reflection 

1.  Terrain,  obstacles.-  The  nature  of  the  terrain  or  ob- 
stacles affects  the  resulting  disturbance.  If  the  terrain  is  smooth,  the 
reflection  appears  to  originate  from  a localized  area  or  point.  If  the  ter- 
rain is  rough,  there  is  more  scattering  and  the  effect  is  not  easily  localized. 
The  terrain  or  structures  becomes  rough  when  the  dimensions  of  its  variations 


^.D.  Kraus,  Antennas.  New  Yorks  McGraw-Hill  Book  Co.,  p,  252,  1950. 


are  comparable  to  the  wavelength  of  the  signal.  It  Is  more  difficult  to 
eliminate  this  source  of  disturbance  than  an  apparent  point  source. 

Obviously,  vegetation  and  other  natural  growth  which  changes  with 
the  seasons  as  well  as  precipitation  (rain  and  snow)  affect  reflections. 
Again,  the  exact  nature  of  the  effect  is  dependent  on  the  wavelength  of  the 
signals. 

2.  Antennas.  - The  pattern  of  the  illuminating  antenna  determines 
what  portion  of  the  surroundings  is  illuminated  from  which  energy  can  be  re- 
radiated. The  aperture  of  the  receiving  (tes")  antenna  determines  what 
portion  of  this  reradiated  energy  is  intercepted.  Furthermore,  all  but  the 
simplest  antennas  are  designed  to  operate  in  specified  frequency  ranges. 

When  operated  outside  their  prescribed  frequency  ranges,  their  radiating 
patterns  and  receiving  apertures  change.  For  example,  the  effective 
aperture,  Ae,  of  an  antenna  is  defined  as  the  power  W,  in  the  terminating 
impedance  of  that  antenna  to  the  power  density,  P,  of  the  incident  wave: 

A --SL 
e P 

2 

It  follows  from  this  definition  that  the  maximum  effective  aperture  of 
a lossless  short  dipole  whose  length,  i , is  much  less  than  the  wavelength, 

X,  of  the  signal  (*.<<!)  is: 

A = 0.119X  2. 
em 

Recall  that  in  free  space  fX=c  where  f is  the  frequency  of  the  signal  and 
c is  the  velocity  of  light. 

When  antennas  are  operated  outside  the  range  for  which  designed,  it 


2 


J.  D.  Kraus,  Antennas , New  York:  McGraw-  Hill  Book  Co.,  Ch.  3,  1950. 
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is  sometimes  difficult  to  obtain  even  approximate  expressions  for  their 
equivalent  apertures. 

3.  Aircraft  mounted  antennas.  - It  is  desired  in  these  tests  to  ob- 
tain the  characteristics  of  the  antennas  in  their  operating  environment  and 
not  simply  their  free-space  patterns.  Consequently,  it  Is  important  to 
isolate  the  reflections  due  to  test  facilities,  terrain,  and  surrounding 
structures  from  those  due  to  the  aircraft  structures  on  which  the  antennas 
will  be  mounted  and  which  will  form  part  of  the  antenna  system. 

By  free-space  pattern  is  meant  the  antenna  pattern  when  an  antenna  is 
isolated  in  a vacuum  and  all  ponderable  bodies  are  an  infinite  distance  from 
the  given  antenna. 


2.3  Determination  of  Sources  of  Reflections 

1.  When  a measured  antenna  pattern  displays  distinct  peaks  and  nulls 
which  are  not  predicted  by  its  design,  the  cause  is  probably  Isolated 
sources  of  reflections. 


The  approximate  location  of  the  source  of  reflections  can  be  determined 

>»■  | Ah  ' 2 \ / , I 


as  follows 


In  Figure  2,  it  is  assumed  that  the  transmitted  signal  is  a vertically1 
polarized,  plane  wave.  Also,  it  is  assumed  that 


where  d is  the  aperture  of  the  receiving  antenna  and 
X is  the  wavelength  of  the  signal. 


At  the  receiving  site,  suppose  a probe  is  moved  vertically  and  the 


measured  field  is  as  shown  In  Figure  3.  For  purposes  of  analysis,  the 
position  of  maximum  measured  field  is  labeled  B and  the  position  of  minimum 
measured  field  is  labeled  C. 

At  position  B,  the  sum  of  the  reflected  signal,  VR,  plus  the  direct 
signal,  Vp, is  a maximum: 

VB  ' V*  ' VD  + h ' V1  * “> 


Note  that  if  Vn  is  given,  V may  be  determined  by  specifying  the 
D R 

boundary  conditions  at  the  point  of  reflection.  Consequently,  VR=  aV^  where 
a is  a complex  constant.  This  relation  has  been  used  in  the  above  equation. 

At  C,  the  sum  of  the  two  signals  is  a minimum 


V_  - V . *■  Vn  - V - V CL  - a) 

C min  D R D 


Incidently : 


V + V , 
max  min 


a “ 


V - V , 
max  min 

V + V . 
max  min 


(2) 

(3) 


If  it  is  assumed  that  the  path  length  of  the  direct  signal  is  approxi- 
mately the  same  for  the  receiving  probe  at  points  B and  C,  the  path  length 
of  the  reflected  signal  at  point  C muat  be  a/2  longer  than  that  at  point  B. 
The  wavelength  of  the  signal  is  X. 


Therefore,  referring  to  Figure  3: 


Cos  0 ” 


2 (Ah)  * 


C4) 


■ • . •»  , v:  ■ 

■ •'Site- 
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RECEIVER 


DIRECTLY  INCIDENT  SIGNAL 


In  other  words,  by  sighting  at  an  angle  0 from  the  receiving  site, 
the  location  of  the  point  of  the  reflection  can  be  determined. 

One  technique  used  to  eliminate  the  effects  of  a reflection  is  the 
installation  of  reflecting  fences.  This  is  illustrated  in  Figure  4, 

Figure  4a  shows  a receiving  antenna  illuminated  by  the  direct  signal 
from  a transmitting  antenna  plus  a reflected  signal.  In  Figure  4b,  a re- 
flecting fence  of  conducting  material  has  been  Installed  at  the  proper  place 
and  orientation  which  intercepts  and  redirects  the  reflected  signal  so  that 
it  does  not  interfere  at  the  receiving  site. 

The  use  of  reflecting  fences  has  been  the  most  successful  method  of 
eliminating  reflections  at  the  Newport  Test  Range,  In  fact,  to  test  a par- 
ticular antenna  at  a given  frequency  this  technique  is  entirely  satisfactory. 
In  practice,  however,  the  task  becomes  prohibitive  if  the  range  is  used  oyer 
a wide  frequency  range  and  for  a variety  of  antennas  because  too  many  fences 
are  required  to  neutralize  all  the  points  of  reflection.  As  a compromise, 
two  different  ranges  were  used  to  cover  the  present  frequency  range  of  in- 
terest. 


£ 

2.4  Results  of  Corrective  Action 

There  are  four  sites  in  operation  at  the  Newport  Test  Range  as  indi- 
cated in  Figure  1 on  Tanner  Hill  and  Irish  Hill.  Three  ranges  are  available. 
Building  1600  on  Tanner  Hill  is  one  terminus  for  each  of  the  ranges.  In 
addition.  Buildings  1625,  1620,  and  an  antenna  tower  711  ft.  west  of  building 

*A11  data  and  figures  labelled  with  an  asterisk  (*)  were  obtained  by  Rome 
Research  Corporation  (RRC)  personnel.  Details  of  their  efforts  are  re- 
corded in  their  monthly  reports. 
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ILLUMINATING  RECEIVING 
ANTENNA  ANTENNA 


o)INTERFERING 

REFLECTION 


b) INTERFERENCE  ELIMINATED 
BY  REFLECTING  FENCE 


Figure  4.  Elimination  of  Reflections  by  Fences. 


1600  on  Tanner  Hill  each  serve  as  the  other  site  for  separate  ranges. 

1.  Range  #2:  Bldg.  1600,  Tanner  Hill  - Bldg.  1625,  Irish  Hill 

An  elevation  profile  of  the  line-of-sight  (LOS)  terrain  be- 
tween buildings  1600  and  1625  is  shown  in  Figure  5.  In  May,  1976,  re- 
flecting fences  as  indicated  in  Figure  6 were  installed.  The  installation 
consists  of  four  fences,  fifty  feet  long  and  four  feet  high.  The  wire  is 
one-half  inch  square  mesh.  The  fences  follow  the  contour  of  the  terrain  and 
are  inclined  uphill  to  reflect  the  transmitted  signals  skyward.  These 
fences  reduced  reflections  in  the  200  MHz  - 400  MHz  frequency  range  but  the 
study  was  incomplete  and  more  work  is  needed  on  this  problem. 


2.  Section  Range:  Bldg.  1600,  Tanner  Hill-Bldg.  1620  Irish 
Hill.  - A profile  sketch  of  the  terrain  between  buildings  1600  and  1620  is 
shown  in  Figure  7. 

The  field  was  probed  as  a function  of  elevation  in  the  0,5  GHz  to 
4.0  GHz  range  on  Tanner  Hill  at  the  building  1600  site  where  it  was  illu- 
minated by  a parabolic  dish  antenna  at  the  Irish  Kill  site.  A typical  plot 
of  the  field  is  shown  in  Figure  8.  Such  a plot  indicates  reflections. 

The  length  of  vertical  travel  from  top  to  bottom  is  16  feet.  If  the 
received  signal  were  a uniform  plane  wave,  the  plot  would  be  a constant. 

The  presence  of  approximately  periodic  peaks  and  nulls  suggests  that  re- 
flections are  present  as  indicated  in  Figure  3.  The  fact  that  the  ratio 
of  peaks  to  nulls  is  not  constant  over  the  entire  plot  suggests  that  more 
than  one  reflection  is  present  and  that  they  do  not  all  contribute  signi- 
ficantly to  the  received  signal  at  all  positions  of  the  probe. 

Further  investigations  revealed  that  when  using  available  antennas. 
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Elevation  Profile  of  Newport  Test  Range  #2 


FENCEPOSTS 


Figure  6.  Layout  of  Reflecting  Fences  at  Test  Site  #2.  Fences  consist  of 
one-half  inch  mesh  in  50  ft.  sections,  four  ft.  high. 

Approximate  locations  of  the  fences  are  shown  in  Figure  22. 


Figure  7:  Profile  Sketch  and  Section  Range. 


Figure  8.  Probe  of  Relative  Power  at  Receiving  Site  of  the  Section  Range  (*Measurements  by  RRC) 


the  total  area  of  the  terrain  Illuminated  was  too  great  to  cover  with  re- 
flecting fences  when  the  Section  Range  was  used  for  signals  covering  the 
frequency  range  of  0.5  GHz  to  4.0  GHz.  Consequently,  It  was  decided  to 
use  this  range  only  for  2.0  GHz  to  4.0  GHz  at  present. 

After  further  study,  reflecting  fences  were  installed  as  shown  in 
Figure  9.  A copy  of  two  plots  in  Figures  10  and  11  show  the  reductions  In 
reflections.  Figure  10  is  a copy  of  the  vertical  pattern  recorded  before 
any  fences  were  erected.  Notice  that  one  sldelobe  is  only  8 db  below  the 
main  beam  and  the  plot  is  highly  unsymmetrical.  Figure  11  shows  the  plot 
after  the  reflecting  fences  were  installed.  Here  the  pattern  more  closely 
approximates  the  free  space  pattern  with  successive  sidelobes  decreasing 
as  predicted. 

3.  Short  Range  Corrections.  - In  checking  the  Section  Range  (Bldg 
1600-Bldg  1620)  it  was  found  that  the  available  antennas  would  not  perform 
satisfactorily  over  the  desired  frequency  range  of  0.5  GHz  to  4.0  GHz.  Re- 
flecting fences  were  installed  which  gave  satisfactory  results  in  the  range 
of  2.0  GHz  to  4.0  GHz  as  described  in  the  above  section.  However,  below 
2.0  GHz  the  beam  width  of  the  antennas  was  so  great  that  a large  area  of 
terrain  was  illuminated  which  resulted  in  interfering  signals  at  the  an- 
tenna site  on  Bldg.  1600.  Installing  additional  fences  was  too  big  an 
effort.  Changing  antennas  during  a test  was  too  unwleldly  and  time  con- 
suming. Consequently,  it  was  decided  to  use  the  so-called  Short  Range  for 
measurements  in  the  0.5  GHz  to  2.0  GHz  frequency  range. 

A sketch  of  this  range  is  shown  in  Figure  12a  and  Figure  12b. 

Figures  13  and  14  show  a probe  of  the  field  (hopefully  uniform)  taken 
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In  initial  installation,  fences  were  four  feet  high. 
In  final  installation,  fences  were  seven  and  one- 
half  feet  high  (eight-feet  fencing  at  twenty  degree 
angle  of  repose) . 

Fencing  consists  of  one-quarter  inch  mesh  wire. 
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Figure  9.  Initial  Reflecting  Fence  Installation  on  Section  Range 
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Figure  12b.  Range  Elevation  Profile,  Short  Range  - Newport  Site 
(*  Measurement  by  RRC) . 


21 
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(TO  RECEIVER  TOWER) 

Figure  15.  Fence  Configuration,  Short  Rang**  Each  fence  is  8 feet  high 
and  sloped  at  approximately  m the  vertical. 

The  fences  consist  of  1/4"  assn  v 
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REMARKS-  2.0  -0.5  GHz 


Probe  of  Relative  Power  at  Receiving  Site  of  Short  Range  with  Reflecting  Fences 
Installed.  (*Data  taken  by  RRC) . 


with  no  modifications  to  eliminate  reflections.  The  plot  should  be  a 
horizontal  line  or  a constant.  Figure  15  shows  a plan  of  the  reflecting 
fences  which  were  added.  Figures  16a  and  b show  the  fields  after  the 
fences  were  in  place. 

The  use  of  fences  here  gave  the  most  dramatic  improvement.  In  Figures 
13  and  14  the  original  field  plots  vary  by  10  and  12  db,  respectively.  After 
installing  the  fences,  when  both  horizontal  and  vertically  polarized  uniform 
plane  wave  signals  were  transmitted,  the  peak-to-peak  variations  of  the  re- 
ceived signals  were  less  than  2 db  over  the  frequency  range  of  0.5  GHz  to 
2.0  GHz.  This  is  illustrated  by  the  plots  in  Figures  16  a and  b. 

The  approximate  location  of  the  installed  fences  is  shown  in  Figure  22 
at  the  end  of  the  report.  This  may  be  compared  with  Figure  1,  the  topo- 
graphic view  of  the  test  site. 

2.5  Cepstrum  Analysis  and  Correction* 

Installing  reflecting  fences  and  absorbing  material  may  be  a major 
effort.  Under  proper  circumstances,  the  true  antenna  pattern  can  be  ex- 
tracted from  the  pattern  distorted  by  reflections  if  the  correct  data  pro- 
cessing is  performed.  This  may  eliminate  the  need  for  installation  of 
fences  or  absorbing  material.  In  addition,  sometimes  after  a major  reflec- 
tion is  eliminated  by  the  installation  of  fences,  lower  level  reflections 
may  be  revealed.  These  additional  effects  can  be  eliminated  by  analytical 
techniques  in  the  data  processing. 


* 

Further  discussion  of  Cepstrum  Analysis  is  Included  in  Appendix 
II,  page  28. 
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The  basis  for  this  analytical  technique  is  the  cepstrum  analysis 
which  will  now  be  developed. 

In  Figure  17,  let 

x(0)  = measured  radiation  pattern  due  to  the  direct 

signal  along  path  R only  (5) 

oix(0-0q)  " measured  pattern  due  to  the  multipath  signal 

along  (R^  + R2).  (6) 

a is  a complex  reflection  factor  and  0^  is  the  angle  between  direc- 
tions R and  R£. 

The  actual  measured  pattern  is  the  sum  of  the  two  above  signals  or 

y (0)  = x(0)  + ox(0  - 0Q)  (7) 

Usually,  | a | <1  and  0^  > 0. 

If  a and  0^  can  be  determined,  and  if  y(0)  is  measured,  Eq.  (7)  can 
be  solved  to  yield  x(0),  the  true  antenna  pattern.  For  simplicity,  only  one 
reflection  is  considered  but  the  technique  can  be  extended  to  include  several 
reflections.  The  solution  of  Equation  (7)  can  be  implemented  also  by  hard- 
ware as  indicated  in  Figure  18.  Note  that  the  input  is  the  total  received 
signal  and  has  been  designated  y(0).  The  output  is  the  desired  antenna  pat- 
tern and  has  been  designated  x(0).  By  inspection  of  the 'circuit  diagram  in 
Figure  18: 

x(0)  - y (0)  - ox(0-0Q). 

The  quantities  a and  0Q  can  be  determined  by  the  techniques  already  in 
use  to  determine  the  locations  of  the  sources  of  reflections. 


*Bogert , Healy,  Tukey,  "The  Quefrency  Analysis  of  Time  Series  for  Echoes," 
Proc.  Sym,  on  Time  Series  Analysis,  M.  Rosenblatt  (ed.),  J.  Wiley,  1963. 


Figure  17. 


X(0) 


z“fl o REPRESENTS  A 
DELAY  OF  Q0  DEGREES 


Figure  18.  Implementation  of  y(6)  * x(0)  ♦ ox(6-0o) 
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However,  they  also  can  be  determined  by  the  method  of  cepstrum  analy- 
sis. Potentially,  this  will  result  In  an  enormous  savings  In  time  and 

I 

effort.  Furthermore,  the  circuit  of  Figure  18  may  be  replaced  by  a direct 
computation.  In  Equation  (7)  If  a and  6q  are  known  and  yC0)  la  measured, 
x(6)  can  be  determined  by 

x(6)  - y(8)  - ax(6  - 0Q)  (8) 

This  Is  developed  in  what  follows. 

2.5.1  Cepstrum  analysis-  In  this  analysis  It  is  assumed 

that  the  pattern  x(0)  of  the  test  antenna  is  given  and  that  the  Illuminating 
source  is  a uniform  plane  wave.  If  there  is  a single  reflection,  the 
measured  signal  Is 


y(8)  - x(0)  + ox(0  - 0Q); 


M < i,  eQ  > o. 


Define  the  Fourier  transform  of  x(0)  as 


w 

Fx(o>)  - | x(0)e‘JW0d0  - F(x); 


b)  - 2tff. 


L^Coj)  - £n(F(u>))  which  Is  the  logarithm  of  F(w). 
Then  the  following  processing  is  performed  on  y(0). 


Fy(w)  - | lx(0)  + oxCfi  - 0oUe-Ja,ed0 


- [1  + ote  jF  (u>) 


1^(0))  - £n(Fy(u))  - ln(l  + ae  °)  + £n(Fx(aj)). 
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y(n)»  x (n ) -!  a*tn-n0);  0(n-l)-0(n)«  8 '•  a s 2 . nos  2 
x ( n ) - y (n)  - 1/2  x(n-2) 

TABLE  3 \ 

y(n)  x (n)  \ Q 

3.06  3.06  \ ^ 

5.66  5.66  \ \ 

8.92  7.39  \ \ 

10.83  8.00  \n*2  \n«l 

11.09  7.39  \ 

9.66  5.66  \ \ 

6.76  3.06  \ 

2.83  0.00  ns3  \ 

1.53  0.00  ^>*-^§^92  ^ \ 


'8.92 
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t3.06 


'10.83 


3.06// 


11.09 
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x(n-na) 
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Figure  19.  Application  of  Cepstrum  Analyais  to/  Antenna  Pattern  Determination 
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If  Ja|  « 1, 


-jU)0 

L (to)  ■»  ae  0 + £n(F  (to) ) 

y x 


Finally,  taking  the  inverse  Fourier  transform  of  (15)  yields: 

+00 


Ly(aj)e^W®  dw 


(15) 


(16) 


+oo 


2tr 


[ae  U + toFx(u))]e:,U)0  dw 


(17) 


or 


a6(0'  - 0 ) + F 1L  (to)  ; <S  (6)=Dirac  Delta  (18) 
x Function 

=1,  0=0 

=0,  otherwise 


a6(0'  - 0Q)  = F'1(Ly(co))  - F'1(Lx(w)) 


(19) 


The  quantity  F ^(Ly(u)))  determined  by  equation  (16)  is  called  the 
cepstrum  of  y(0).  It  is  the  Inverse  transform  of  the  logarithm  of  the 


direct  Fourier  transform  of  y(0).  Because  of  the  step  of  taking  the  loga- 
rithm of  Fy((j),  the  domain  of  the  independent  variable  of  Equation  (16)  is 
not  9;  it  has  been  named  the  cepstrum.  However,  because  of  the  general  re- 
lationship log  AB  “ log  A + log  B,  the  value  of  a and  0^  can  be  obtained 
from  the  cepstrum  of  a multipath  signal  (or  function)  if  the  undistorted 
value  of  the  signal  is  known. 

Notice  that  this  Inverse  transform  includes  an  impulse  (or  Delta 
Function)  which  is  delayed  0^  and  of  amplitude  a which  allows  determination 
of  these  two  quantities. 

Repeating,  knowing  a and  0^,  the  pattern  x(0)  of  an  unknown  antenna 
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can  be  obtained  by  computation  using  Equation  (8)  or  processing  with  the 
hardware  of  Figure  18. 

Equation  (19)  reveals  why  we  require  the  field  strength  plot  of  an 

antenna  whose  free  space  pattern  Is  known.  Observe  that  magnitude  and 

delay  of  the  reflection  equals  the  cep strum  of  the  received  data  minus 

the  cepstrum  of  the  free  space  or  undlstorted  pattern. 

Equation  (15)  assumes  « 1.  If  this  Is  not  true,  other  Impulses 

will  appear  in  the  solution,  but  they  will  have  values  of  — *■—  6(0  - n0A) 

n!  u 

and  can  be  recognized.  Impulses  due  to  other  reflections  which  appear 
will  not  be  so  related.  Again,  6(0)-l,  0-0  and  6(0)-O,  otherwise. 

Finally,  Equation  (8)  Is  a first-order  difference  equation  and  to 
obtain  a unique  solution,  one  Initial  value  Is  required. 

An  example  of  the  application  of  the  technique  follows  In  the  next 
section. 


2.5.2  Cepstrum  Processing  - The  application  of  cepstrum  processing 
can  be  Illustrated  by  an  example. 

Let  the  solid  line  In  Figure  19  represent  a recorded  field  strength 
plot.  Furthermore,  suppose  that  we  observe  the  data  y(n)  at  discrete 
points  as  a function  of  0,  every  ir/8  radians,  letting  the  positive  vertical 
axis  correspond  to  9 - 0(n  - 0),  and  the  counterclockwise  direction  as 
positive. 

The  observed  values  of  y(n)  are  shown  below  In  Table  1, 


Table  1:  Distorted  Pattern  Values 


0 1 


3 


8 9 10-15 


y (n)  0 3.061  5.657  8.922  10.828  11.087  9.657  6. "57  2.828  1.531  0 


In  addition,  suppose  that  the  desired  free  space  pattern  is  known 


to  be: 


p = -8  sin  0 
= 0 


, 0 0 <_  IT 

, otherwise. 


as  shown  by  the  dotted  circle  in  Figure  19.  Its  pattern  values  are: 


Table  2:  Free  Space  Pattern  Values 


n 

0 

1 

2 

3 

4 

5 6 

7 

8-15 

x(n) 

0 

3.061 

5.657 

7.391 

8 

7,391  5.657 

3.061 

O 

Using  the  data  in  Tables  1 and  2,  the  constants  a and  n^  used  in 
equation  20  will  be  evaluated. 

y (n)  - x(n.)  + ax(n  - nQ).  (20) 

With  a and  n^,  the  effects  of  the  reflection  can  be  eliminated  in 
the  measurements  of  unknown  antenna  patterns. 

These  quantities  will  now  be  evaluated.  Because  we  deal  with 
discrete  data,  the  discrete  Fourier  transform  will  be  used.*  The  cepstrum 
of  x(n)  can  be  computed. 

Oppenheim  and  Schafer,  Digital  Signal  Processing.  Englewood  Cliffs,  N.J.: 

Prentice  Hall,  Ch.  3. 
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XGO  - (Y)x(n)G"jC2lT/N)nK  ; n - 0,1,2, ,..15;  N - 16. 
n-0 

in  X(K)  - in(.  l x(n)e“3  ^8^nK);  n - 0,1, 2, 3 ...  15 
n*0 


x»(n)  - l dn(X(K))e+j(7r/8)Kn;  K - 0,1,2 15 

K-0 


= cepstrum  of  x(n). 


Next  calculate  the  cepstrum  of  y(h), 

(N-D  _w. 


Y (K)  » l y(n)e"3(2Tr/N)nK;  n - 0,1, 2, 3 ...15;  N - 16 
n-0 


- I [x(n)  + ox(n  - n0)]e~jCll/8)nK 
n-0 

-j(2ir/N)nnk 

- X(K) (1  + ae  °) 

-j(2rr/N;n0k 

Jln(Y(K))  - Hn(X(K))  + Jln(l  + ae  ) 

-j  (2ir/N)n  k 

* f,n(X(K))  + ae  ; |a|  « 1. 


Therefore 


y'(n)  - l An(Y'(K)e3<1r/8)Kn 
k-0 


- x' (n)  + a6(n  - nQ) 


a6(n  - nQ)  - y'(n)  - x'(n) 


In  this  particular  case,  calculations  yield: 


a - 0.5, 


(22a) 


nQ  - 2. 


(22b) 


« • - 
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In  other  words,  given  an  antenna  whose  free-space  pattern  x(n)  is  known, 
measure  its  response  y(n)  on  the  range  where  a reflection  may  exist.  Using 
the  calculations  outlined  above,  determine  the  constants  a and  nQ.  The 
results  labeled  (22a),  (22b)  indicate  that  they  will  appear  as  a data  point 
or  impulse  of  magnitude  a at  point  n^.  If  more  than  one  reflection  exists, 
additional  impulses  will  appear  in  (21). 

In  practice,  these  calculations  will  be  done  on  a computer.  However, 
a sample  calculation  can  be  accomplished  with  a hand  calculator. 

The  use  of  this  technique  is  subject  to  the  same  limitation  as  fence 
installations  and  absorbing  material.  Changes  in  frequency  range,  etc. 
mav  introduce  additional  reflections  and  require  a new  determination  of  the 
a's  and  n^'s. 

With  these  constants  for  a given  antenna  range,  the  true  free-space 
pattern  of  an  antenna  can  be  calculated  from  its  measured  pattern  and  the 
use  of  (20),  rearranged. 

Again,  assume  that  the  pattern  shown  as  the  solid  line  in  Figure  19 
has  been  measured.  Furthermore,  assume  that  the  range  has  been  calibrated 
and  that 

a ■ 0.5, 

n0  * 2* 

These  quantities  can  be  determined  by  the  geometric  optical  method 
or  the  cepstrum  method  as  discussed  above. 

With  a and  n^  known,  x(n)  can  be  determined  from  the  above  equation; 


I 


I 


* 


I 


! 
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NORMALIZED  SIGNAL  1 

l 


NORMALIZED  SIGNAL  1 


i.e. , 

x(n)  - y(n)  - ax(n  - nQ) 

* y(n)  - 0.5x(n  - 2). 

Application  of  this  equation  yields  the  value  of  x(n),  the  true 
pattern,  shown  in  Table  3 in  Figure  19.  This  pattern  is  shown  as  the 
dotted-line  circle,  p - -2r  sin  9 - -8  sin  9,  in  Figure  19.  These  are 

the  same  data  as  listed  in  Table  2,  page  14. 

2.6  Alternate  Technique  of  Signal  Processing 

The  example  in  the  last  section  illustrates  a method  of  resolving  a 
scalar  multipath  signal  into  its  direct  and  reflected  components.  This 
method  may  not  be  applicable  to  vector  signals.  However,  some  classes  of 
vector  signals  can  be  analyzed  in  an  analogous  way. 

The  example  below  demonstrates  a technique  of  determining  the  source 
(location  and  magnitude)  of  a disturbing  reflection  using  the  Fourier  trans- 
form alone.  Figure  20  describes  the  situation.  The  electric  field  is 
probed  in  the  vertical  (Z)  direction  and  an  approximately  periodic  pattern 
such  as  Figure  8 is  recorded. 

Assume  that  at  the  origin  the  received  signal  is  a maximum,  i.e.,  the 
direct  and  reflected  signal  are  in  phase  giving 

y(o)  - 1 + r. 

At  any  position  n,  the  path  length  of  the  directly  Incident  signal  is 
approximately  the  same  but  the  reflected  signal  has  traveled  an  added  dis- 
tance: 
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Ad  - (n)  (sin  8) . 

The  received  signal  of  wavelength  A may  be  written  as: 

-3(Z  /A)  (27t/N)  (n)  (sin8) 


y (n)  - 1 + re 


max 


1 + re 


-jk  (2ir/N)n 


k - sin  8. 

o A 


(23) 


Recall  that  the  discrete  Fourier  transform  is  defined  as: 

(N-l)  -j  (2ir/N)kn 
Y(k)  - l y(n)e 
n-0 


: J 


(N-l)  -j  (2ir/N)  (k_)  (n)  -j (2u/N)kn 

_ I l (1  ♦ re  u )e 

" N n-0 

- 6(k)  + r6(k  + kQ)  * 6(k)-l,  k-0  and  6(k)-0  otherwise.  (24) 

Examination  of  this  result  reveals  that  the  Fourier  transform  of  y(n) 
consists  of  two  impulses.  The  impulse  at  the  origin  is  due  to  the  direct 
signal  and  the  delayed  Impulse  results  from  the  reflection.  In  this 
example  the  direct  signal  has  been  assumed  to  have  an  amplitude  of  one. 

The  ratio  of  the  magnitudes  of  the  two  impulses  gives  the  relative  vslue  of 
the  reflection.  This  quantity  has  been  denoted  by  a in  the  previous  problem. 

The  delay,  kQ,  allows  determination  of  the  direction  (and,  in  turn,  the 
location)  of  the  reflection  from  the  earlier  formula 

k -is  sin  8 

O A 


or 


k A 


sin  8 - 


max 
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In  fact,  when  r and  kQ  have  been  determined,  the  true  antenna  pattern 
can  be  determined. 

Suppose  r and  kQ  are  known,  we  may  represent  an  arbitrary  antenna 
pattern  as: 

-j(2TT/N)k  n 
y(n)  - x(n) [1  + re  ] 

where  x(n)  is  the  pattern  if  no  reflection  were  present. 

The  magnitude  of  the  measurement  may  be  represented  as: 


| y (n)  | = y x^(n)f(l+T  cos  ~ nk  )2  + r sin  nk 

» NO  NO*4 

All  the  quantities  in  this  expression  are  known  except  x(n)  which  can 
be  evaluated  at  every  data  point,  n. 

To  illustrate  this  method  it  is  applied  to  the  data  recorded  in  Fig.  8. 
Note  that  the  units  of  this  data  are  decibels  (db).  An  arbitrary  reference 
level  was  assumed  and  the  data  was  converted  to  the  units  of  voltage. 
Following  this,  the  Fast-Fourier  Transform  was  taken.  The  result  is  shown 
in  the  table  below. 

Examining  Y(k)  in  Table  4 shows  that  there  is  reflection  corresponding 

to  k - k ■ 5,  The  value  of  the  impulse  at  k-5  is  1.46.  However,  the  value 
o 

of  the  Impulse  at  k*0  in  Table  4 due  to  the  direct  signal  is  8.92. 

Therefore,  the  relative  magnitude  of  the  reflection,  rj,  is  the  ratio  of 
these  two  numbers.  This  is  a specific  value  for  the  quantity  r shown,  in 
equation  (23)  and  (24). 


Its  relative  magnitude  is: 

1.46  A 

rl  " 8.92  0,16 

Its  direction  is  determined  by  the  formula: 


■ 

* 

) 


< 
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Table  4.  Electric  Field  Probe  Data  and  Its  FFT 


(from  Figure  8) 


n y (n)  | Y(k)  | - |±  (n)e~j  (27r/N)kn | 

n-0 


0 

5.82 

8.92 

0 

1 

8.91 

0.32 

1 

2 

12.02 

0.17 

2 

3 

9.33 

0.17 

3 

4 

6.10 

0.82 

4 

5 

11.89 

1.46 

5 

6 

11.89 

0.23 

6 

7 

3.98 

0.31 

7 

8 

10.0 

0.11 

8 

9 

13.03 

0.31 

9 

10 

7.08 

0.23 

10 

11 

7.08 

1.46 

11 

12 

10.0 

0.82 

12 

13 

7.94 

0.17 

13 

14 

7.59 

0.17 

14 

15 

10.12 

0.32 

15 

k - sin  8 


Ak 


sin  8 - - — - ^ * 1 » Z - 16'. 

Z 16  max 

max 


8 - 18. 2e 
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Using  the  optical  method  described  earlier  in  the  report.  It  was 
determined  that  the  major  reflections  came  from  a direction  bracketed  by 


.*S**#v- 


angles  of  16.7°  to  18.1  degrees. 


Further  examination  of  Table  4 indicates  that  there  is  a reflection 
corresponding  to  k - 4.  Its  magnitude  is  ■ 0.09.  The  recorded  data  in 
Figure  8 indicate  that  there  is  more  than  one  reflection.  In  order  to  in- 
terpret the  remaining  points  in  Y(k),  further  study  is  required.  However, 
the  technique  does  show  enough  promise  to  warrant  that  study. 

2.7  Conclusion 

Elimination  of  the  effects  of  reflecting  signals  is  necessary  on  out- 
door antenna  ranges  if  meaningful  pattern  data  is  to  be  obtained. 

The  proper  use  of  reflecting  fences  can  be  very  effective  in  eliminat- 
ing reflections.  However,  a particular  installation  may  be  useful  for  only 
a limited  frequency  range.  This  is  because  usually  the  beam  width  or  aper- 
ture of  the  illuminating  antenna  varies  with  frequency,  and  thus  exposes 
more  (or  less)  of  the  surroundings  to  radiation,  giving  rise  to  additional 
reflections . 

Providing  fences  to  be  effective  for  the  general  case  would  require 
an  effort  which  may  be  practically  prohibitive.  If  the  beam  width  or  pat- 
tern of  the  illuminating  antenna  changes  (say,  because  of  mechanical  mis- 
alignment of  its  elements  or  changes  in  frequency  of  the  input  signal), 
additional  reflections  will  probably  be  introduced.  Therefore,  it  is  im- 
portant to  invest  in  illuminating  antennas  which  meet  the  specifications 
for  which  the  constructed  fences  were  designed  and  insure  that  they  are 
installed  properly. 

In  addition,  absorbing  material  also  can  eliminate  reflections.  It  is 
most  useful  in  eliminating  reflections  from  portions  of  the  receiving 

I 

structure.  It  is  too  expensive  to  replace  fences  eliminating  ground  re- 
flections and  offers  no  particular  advantages. 
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Finally,  the  elimination  of  the  effects  of  reflections  by  processing 
the  observed  data  would  represent  a large  saving  in  time  and  effort.  Most 

of  the  equipment  to  Implement  this  technique  is  already  in  place.  Data 
processing  equipment  is  already  used  to  process  the  measurements.  Also, 
the  range  equipment  is  calibrated  every  day  before  measurements  begin, 
using  standard  gain  horns.  These  could  be  used  for  cepstrum  determina- 
tion. 

Furthermore,  the  particular  method  outlined  here  is  only  one  of 
several  based  on  this  general  technique.  For  Instance,  the  antenna  could 
be  fixed  in  place  and  a known  signal  as  a function  of  time  recorded  at  the 
receiver.  This  time  dependent  data  could  be  processed  rather  than  the 
spatial  data  to  yield  the  same  information.  Practical  considerations  would 
decide  which  specific  method  to  Implement.  The  application  of  cepstrum 
analysis  to  this  problem  as  described  in  this  report  shows  promise  and 
should  be  pursued. 


' 

' I 
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3.  INCREASED  FREQUENCY  COVERAGE:  18  GHz  - 40  GHz 


3.1  Introduction 

The  IEEE,  in  its  recent  Standard  Letter  Designations  for  Radar  Frequency 
Bands,  has  designated  18-27  GHz  as  K band  and  27-40  GHz  as  Ka  band.  Note  that 
these  frequencies  correspond  to  free-space  wavelengths  of  1.67  cm  to  0.75  cm. 
One-quarter  wavelength  of  a 40  GHz  signal  is  less  than  2mm  in  length.  Vibra- 
tions or  mechanical  variations  of  an  aircraft  structure  of  this  order  of  mag- 
nitude can  have  profound  effects  on  the  performance  of  antenna  mounted  on  the 
aircraft  unless  they  are  designed  to  be  immune  to  these  effects. 

Also,  there  are  well-known  oxygen  and  water  absorption  lines  above  20 
GHz  . This  causes  atmospheric  transmission  characteristics  to  vary  with 
weather.  These  data  should  be  obtained  before  beginning  antenna  tests  in 
this  frequency  range.  See  Figure  21. 

3.2  Equipment : Receivers 

The  Scientific-Atlanta,  Inc.  series  1711  Microwave  Measurement  Receiver 
covers  the  frequency  range  of  2.0  GHz  to  40.0  GHz.  The  receiver  has  a dyna- 
mic range  of  60  db  and  a sensitivity  of  (-75)  dbm  in  the  20-40  GHz  frequency 
range.  A series  of  RF  to  IF  waveguide  mixers  is  used  to  cover  this  range. 
Model  (13A-18)  Tunable  Waveguide  mixer  covers  18-26.5  GHz  and  model  (13A-26) 
covers  26. 5 -49.0GHz. 

In  fact,  the  frequency  range  of  the  receivers  in  use  at  the  Newport 
sites  can  be  extended  to  40  GHz  by  the  addition  of  these  waveguide  mixers  . 
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3.3  Equipment:  Signal  Sources 

Several  manufacturers  market  signal  generators  in  this  frequency  range 

* 

with  power  outputs  up  to  the  order  of  100  mw.  Higher  power  sources  do  not 
seem  to  be  available  (1977). 

However,  traveling-wave  tube  amplifiers  are  available  with  power  outputs 
of  over  1 kW  (e.g.  from  Hughes  Aircraft  Company).  There  are  a number  of 
radar  systems  operating  in  this  frequency  range  and  there  are  many  device, 
component,  and  system  manufacturers  active  in  the  millimeter  wave  area.^ 

Listed  below  are  the  specifications  of  two  instruments  currently  avail- 
able. 


Manufacturer 


Watkins- Johnson  Co. 


Hewlett  Packard  Co. 


Model  Number  Frequency  Range  Power  Output 


WJ-1250 


938A 

940A 


18-26.5  GHz 
26.5-40  GHz 

18-26.5  GHz 
26.5-40  GHz 


10  mW 
5 mW 

Output  depends  on 
input.  0. 5-1.0  mW 
(typical) 


Note:  The  H-P  938A  and  940A  are  frequency  doublers  and  require 
input  drivers. 

The  required  power  output  of  the  signal  source  can  be  estimated  as 
follows. 

P - S + N - Gr  + tJ  - Gt 

where  P ■ Required  power  output  of  signal  source  in  dbm. 

S - Sensitivity  of  receiver  in  dbm. 

N - Variations  of  maxima  to  nulls  in  pattern  in  db. 


4. 

See  Microwave  Journal,  Vol.  20,  No.  11,  November  1977 

See  Reference  Data  for  Radio  Engineers.  4th  edition,  ITT  Corp.  Fig.  32, 
page  752.  Fig.  21  of  this  report  indicates  that  loss  due  to  water  vapor  and 
oxygen  absorption  will  be  less  than  1 db. 
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G ■ Receiver  antenna  gain  in  db. 

K 

T * Transmission  attenuation  in  db. 

L 

G„  = Transmitter  antenna  gain  in  db. 

T 

Some  estimates  of  required  power  are  listed  in  the  table  below. 

The  smallest  estimate  is  obtained  if  a frequency  of  20  GHz  is  considered 

and  the  Short  Range  (distance  about  0.25  km.)  is  used.  In  all  calculations, 

S = -75  dbm  and  G = 0 db.  In  this  particular  case  T = 109  db;  let  N = 

R u 

30  db  and  GT  = 30  db.  Then: 

P = S + N - Gjj  + T^  - G^.  ; (Short  Range,  20  GHz) 

= -75  + 30  + 0 + 109  - 30 
= 34  dbm 


= 2.5  watts 

The  largest  power  estimate  is  obtained  if  40  GHz  signal  sources  on  Range 
#2  (distance  about  2.5  km)  with  N = 40  db  and  GT  * 30  db  are  considered. 

Then:  P = S + N - GR  + T^  - G^  ; (Range  //2,  40  GHz) 

- -75  + 40  + 0 + 132  - 30 

= 67  : 70  dbm 
L 

- (10)  watts 

The  gain  of  the  transmitting  antenna  is  a parameter  under  the  designer's  con- 
trol and, obviously,  is  very  important.  Here  GT  is  assumed  to  be  30  db.  Any 
change  in  its  value  directly  affects  the  required  power  of  the  signal  source. 

The  published  sensitivity  specification  in  the  26  to  40  GHz  range  for 
the  Scientiflc-Atlanta,  Inc.  receiver,  series  1740,  which  employs  phase-lock 
techniques,  is-100  dbm.  Of  course,  using  this  receiver  sensitivity  in  the 
above  calculations  reduces  the  required  signal  power  by  25  db.  However,  in 
addition,  a phase-lock  reference  input  that  is  35  db  greater  than  the  signal 
channel  is  needed. 
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Table  5 Required  Signal  Power 


Range 

Freq. 

S 

N 

gr 

tl 

gt 

P (dbm) 

P (watts) 

Short 

20 

75 

30 

0 

109 

30 

34 

2.5 

Short 

40 

75 

40 

0 

115 

30 

50 

100 

n 

20 

75 

30 

0 

126 

30 

50 

100 

n 

40 

75 

40 

0 

132 

30 

90 

(10)4 

3.4  Conclusions 

The  short  wave  lengths  corresponding  to  this  frequency  range  require 
Increased  attention  to  mechanical  details  in  making  measurements. 

The  absorption  by  oxygen  and  water  of  electromagnetic  energy  in  this 
frequency  range  causes  increased  atmospheric  transmission  loss. 

Receiving  equipment  in  this  frequency  range  is  available.  Transmitting 
equipment  may  not  be  available  off  the  shelf,  but  it  can  be  provided  on 
special  order. 
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Attenuation  in  Decibels  per  Kilometer 


Figure  21.  Atmospheric  Absorption  by  the  1.35  cm  Line  of  Water  Vapor  and  the 
0.5  cm  Line  of  Oxygen.  Water  Vapor  Absorption  Calculated  for  the 
Mean  Absolute  Humidity  at  Washington,  D.C.  (From  NBS  Report 
6763,  Attenuation  of  Microwaves,  by  B.R.  Bean  and  L.P.  Riggs, 
April  21,  1961) 


IRISH  HILL 
BLDG.  1625 
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REPORTS  OF  EARLIER  RADC  STUDIES 


1.  "Investigation  of  Precision  Antenna  Pattern  Recording  and  Display  Tech- 
niques," Final  Report,  AD415-912,  March,  1963. 

2.  "Investigation  of  Precision  Antenna  Pattern  Recording  and  Display 
Techniques,  Phase  II,  Final  Report,  J.  Searcy  Hollis  et  al,  AD630-124, 
February  1966. 

3.  "Error  Analysis  for  the  RADC  F-4  Antenna  Test  Facility,  Final  Report, 
Mario  Macera,  RADC-TR-74-309,  December  1974,  B001090L. 
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Appendix  II 

An  Example  of  Cepstrum  Analysis 

> Elimination  of  Effects  of  Reflections  by  Data  Processing 

I.  Introduction 

Suppose  that  the  signal  y(n)  shown  In  Fig.  A is  received  or  measured. 
Further,  suppose  that,  in  fact,  y(n)  Is  really  the  sum  of  two  signals. 


x(n)  + w(n),  shown  in  Figs.  B and  C. 

Anally,  suppose  that  w(n)  is  an  attenuated  and  delayed  replica  of 
x(n);  i.e., 

w(n)  - a x(n-nQ) 

- 0.4  x(n-2) 

f 

Then:  y(n)  ■ x(n)  + a x(n-n^) 

■ x(n)  + 0.4  x(n-2) 

Therefore  x(n)  • y(n)  - 0.4  x(n-2). 

With  the  y(n)'s  known,  if  a single  value  of  x(n)  Is  known,  all  the  ocher 
values  of  x(n)  can  be  determined  by  the  above  formula. 

Suppose  it  Is  known  that  y(0)  - x(0)  - 1,  (again,  as  indicated  in 
Figs.  A and  B).  Then  it  follows: 


J.-.r  * V>^'- -V.~  V*.*  ' 


n y(n) 

0 1.0 

1 0.5 

2 1.4 

3 0.7 

4 1.4 

5 0.2 

6 0.4 

n > 6 0 


x(n)  - y(n)  - 0.4  x(n-2) 
1.0  -1.0-0 
0.5  - 0.5  - 0 
1.0  - 1.4  - 0. 4(1.0) 

0.5  - 0.7  - 0.4(0. 5) 

1.0  - 1.4  - 0. 4(1.0) 

0 - 0.2  - 0.4(0. 5) 

0 - 0.4  - 0. 4(1.0) 
0-0-0 


Of  course,  chese  values  of  x(n)  agree  with  those  of  Fig.  B. 

Note:  The  Independent  variable  may  represent  time.  Then  the  received 
signal  y(n)  results  when  a transmitted  signal  x(n)  is  distorted  by  a multi- 
path  transmission  link.  In  another  case,  the  Independent  variable  may 
represent  a physical  dimension  such  as  elevation  or  azimuth.  Then  y(n) 
may  represent  a received  signal  which  is  composed  of  a directly  trans- 
mitted signal  x(n)  plus  w(n)  - a x(n-n^) , its  reflection  from  some  local 
obstacle. 

The  development  is  presented  for  discrete  data  but  applies  equally  as 
well  to  continuous  functions.  Also,  for  simplicity  the  discussion  con- 
siders only  one  delayed  (or  distorting)  signal;  however,  the  analysis  may 
be  extended  to  include  multiple  delays. 

2.  Determination  of  a and  n by  Cepstrum  Analysis 
o 

The  above  example  was  fabricated  end  since  a and  nQ  were  assumed, 
x(n)  was  easily  recovered  when  y(n)  was  given. 


54 


A natural  question  to  ask  Is  the  following:  Can  a given  trans- 
mission link  or  antenna  range  be  "calibrated"?  Can  a and  n^  be  deter- 
mined so  that  any  unknown  x(n)  can  be  recovered  when  y(n)  is  received? 

Fortunately,  the  answer  to  this  question  is  yes.  If  a known  signal, 
x(n) , is  transmitted  and  the  received  signal,  y(n),  is  recorded,  the  tech- 
nique of  cepstrum  analysis  will  yield  a and  nQ. 

Again,  consider  the  equation: 


y(n)  ■ x(n)  + a x(n-nQ) 


In  this  case,  the  y(n)'s  and  x(n)'s  are  known.  The  constants  a and  nQ 
are  to  be  determined.  This  can  be  done  as  follows: 

First,  take  the  Discrete  Fourier  Transform 

(»-l)  -J  n k 

Y(k)  - l y(n)  e 
n*0 


2ir 


-J  £*k 


<J"1)  , , (H-l)  , „ 

- I x(n)  c N + 0 l x(n^  ) e » 

n“°  n-0  ° 


A— 

- X(k)  ( 1 + a c"3  N “o) 


where  Y(k)  is  the  Discrete  Fourier  Transform  of  the  function  y(n)  and  X(k) 


*Bogert,  Healy,  Tukey,  "The  Quefrency  Analysis-  of  Time  Series  for 
Echoes",  Proc.  Symposium  on  Time  Series  Analysis.,  M.  Rosenblatt  (ed.), 
J.  Wiley,  1963. 
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is  the  Discrete  Fourier  Transform  of  x(n).  N represents  the  total  number  of 


points  of  data 


Next,  take  the  logarithm  of  each  side 


ln(Y(k))  - ln(X(k))+  ln(l  + a e 


ln(X(k))  + a e 


Note:  The  assumption,  a « 1,  is  made  to  simplify  and  clarify  the 


analysis.  If  a is  not  small,  this  development  still  applies 
Finally,  take  the  inverse  Discrete  Fourier  Transform 


1 — 

ln(Y(k))  eJ  N 


[ln(X(k)>  + a c 


The  left-hand  side  of  the  equation  is  the  desired  result.  It  is  an 


impulse  of  strength  a located  at  position  n^.  Because  y(n)  and  x(n)  are 
known,  y'(n)  and  x'(n)  can  be  calculated,  and  a and  n have  been  deter- 


After  a and  n have  been  determined  for  a particular  experiment 


ty  go  back  to  the  technique  illustrated  at  the  beginning  of  this 


IHfetiNiMlISM 


discussion  and  extract  an  unknown  transmitted  signal,  x(n) , from  the 


received  signal  y(n);  i.e., 


x(n)  - y(n)  - a x(n-n  ). 

o 

3.  Discussion 

By  examining  the  steps  in  the  Cepstrum  Analysis,  insight  is  gained 
into  why  it  works. 

We  begin  with: 

y(n)  ■ x(n)  + a x(n-n  ) 

o 

The  right-hand  side  of  this  equation  may  also  be  written  as  the  con- 
volution of  the  signal  x(n)  and  the  characteristic  function  of  the 
transmission  link;  i.e., 

y(n)  - x(n)  * h(n) 

When  the  Fourier  Transform  is  taken,  the  right-hand  becomes  a product;  i.e., 

F(y(n))  - Y(K)  - X(K)  H(K) 

The  Fourier  Transform  is  taken  because  of  the  form  of  H(K) . Recall  if 
the  Fourier  Transform  of  x(n)  is  X(K),  the  Fourier  Transform  of  x(n-nQ)  ■ 
X(K)  « no.  Now  the  inverse  Fourier  Transform  of  ( ^ 8o  ii 

an  impulse  at  n which  can  be  identified  and  will  be  useful. 
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However,  this  quantity  is  submerged  in  the  product  X(K)  H(K) . But 
the  In  XH  ■ In  X + In  H and  the  product  is  transformed  into  a sum;  i.e. 

In  (Y(K))  - In  (X(K))  + In  (H(K))) 

Taking  the  inverse  transform  and  rearranging  it,  the  useful  result  emerges 
as 

F In (Y(k))  » y'(n)  - x'(n)  + a 6(n-n  ) 

o 

a <5(n-nQ)  * y'(n)  - x>(n) 

This  technique  is  illustrated  in  the  example  which  follows.  Con- 
sider a received  signal: 

y(n)  - |x(n)  + ao  x(n-no)  + xfci-n^  | . where  |x|  is  the 

magnitude  of  x. 

- |x(n)  + 0.1  x(n-5)  + j(0.1)  x(n-8) | 

Further,  it  is  known  that: 

x(n)  - 1 , n-0,1,2,3,4 

■ 0 otherwise 

To  be  more  general,  two  reflections  or  delays  are  included,  one  of 
which  has  imaginary  values.  The  data  y(n)  and  x(n)  are  shown  in  Fig.  D 
and  E. 

When  this  data  is  processed  to  ascertain  the  a^'s  and  the  n^'s,  the 
results  are  as  shown  in  Fig.  F.  Note  the  - 0.1  at  nQ-5  and  c^-0.1 


at  1^*8. 

Also,  note  the  unwanted  values  at  points:  2nQ«lO,  nQ  + n^*13, 

2^*16.  These  additional  points  result  from  the  fact  that: 

. ...  x 12.13.  . t , \ ( 1+n ) o£+— ; -1  < o < 1;  n «=  1,2,3— 

ln(l+a)  « b - j#  + a + + _1  ~ 

In  the  discussion,  the  approximation  was  made  that  ln(l+a)  * or,  o<<1. 

Because  these  spurious  responses  occur  at  multiples  of  nQ,  n^  and 
(p  »o  + q n1),  where  p and  q are  integers,  they  can  he  identified  and 
ignored. 
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